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AdS/QCD correspondence predicts the structure of quark-antiquark potential in the static limit. 
We use this piece of information together with the Salpeter equation (Schrodinger equation with 
relativistic kinematics) and a short range hyperfine splitting potential to determine quark masses 
and the quark potential parameters from the meson spectrum. The agreement between theory and 
experimental data is satisfactory provided one considers only mesons comprising at least one heavy 
quark. We use the same potential (in the one-gluon-exchange approximation) and these data to 
estimate constituent diquark masses. Using these results as an input we compute tetraquark masses 
using a diquark-antidiquark model. The masses of the states X(3872) or Y(3940) are predicted 
rather accurately. We also compute tetraquark masses with open charm and strangeness. Our result 
is that tetraquark candidates such as D s (2317), D s (2457) or X(2632) can hardly be interpreted as 
diquark-antidiquark states within the present approach. 

PACS numbers: 12.39.Ki, 12.39.Pn, 12.90.+b 
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I. INTRODUCTION 



Hadron spectroscopy has received in the last few years a renewed attention, both experimental and theoret- 
ical. On the experimental side several new charmonium and charmed states have been observed (for reviews 
see e.g. [l|, Q). The interpretation of the new states has triggered a considerable amount of theoretical work, 
especially so because some of the new states cannot be interpreted as ordinary quark-antiquark mesons (for 
reviews see 0, Q). For some of the new states, such as the A(3872) or F(3940) [E [fi , H B, B , ToL 11 1 an 
interpretation in terms of exotica (diquark-antidiquark bound states) has been given [l2j, |l3|, [lj, [l5|, using a 



refined version of the constituent quark model [16| . On the other hand also more conventional interpretations 
have been proposed. It is therefore useful to have an independent approach to the calculation of tetraquark 
masses in order to assess the validity of the diquark-antidiquark interpretation of the new states. This is 
one of the aims of the present paper. We approach this problem using a semirelativistic method based 
on the use of a relativistic wave equation, the Salpeter equation. This is a Schrodinger equation with the 
relativistic kinematics (for previous use of this equation for hadron spectroscopy see e.g. [17l | and references 
therein); as such, it has the usual limitations of the static potential approach, but presents the advantage of 
the relativistic kinematics. 

A crucial point in this a ppr oach is the choice of the potential. There are several proposals in the literature, 
e.g. the Cornell potential |l8| or its variants 2(j. We shall use here a modified version of a static quark 



potential that has been recently found in the context of the AdS/QCD correspondence 2l|. As shown in 



2l| . this potential has the same behavior expected from QCD, ie it is linearly rising at large distances while 
presenting a Coulomb behavior at small distances. Clearly the interpolation between these two behaviors is 
phcnomcnologically relevant if it corresponds to length scales typical of the hadrons. This is a possible way 
to discriminate among the different models. Another aim of this paper is therefore to find the region where 
the AdS inspired potential differs from the usual QCD-based potentials and to see if it is of phenomcnological 
significance. We fix the parameters of the model using, as an input, data from the meson spectrum. We 
compare our results with the available experimental data for the mesonic spectrum and we find a reasonable 
good agreement. Once the parameters of the model are fixed, we can compute the diquark masses, ie a 
set of phenomcnological parameters that we use, together with some additional hypothesis, to predict other 
spectra. 

The results of the paper are as follows. Assuming a diquark-antidiquark structure of tetraquark states, we 
find masses for some of the X and Y states in reasonable agreement with experiment. The same procedure 
can be applied to tetraquarks comprising one charm quark and hidden strangeness. In this case the results are 



hardly compatible with the masses of possible tetraquark candidates, ie the states D s (2317) 22j, D s (2457) 
[il] and A (2632) 0. 

The paper is organized as follows. In section [TT] we describe the model, introducing the Salpeter equation, 
the potential term and the numerical method adopted to solve the wave equation. In section HTT1 wc deter- 
mine the parameters of our model fitting meson spectra; the results allow a comparison with the presently 
available meson data. In section [TVl the constituent diquark masses are calculated, in the one-gluon exchange 
approximation. In section [V] wc compute the spectra of a few tetraquark states. Finally, in section I VII wc 
draw our conclusions. 
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II. WAVE EQUATION AND THE CENTRAL POTENTIAL 



The wave equation we will use in this paper is an eigenvalue equation for the hamiltonian of two point-like 
particles, taking into account the reltivistic kinematics. It is known in the literature as Salpeter equation. 
In QCD it arises from the Bcthc-Salpetcr equation replacing the interaction with an instantaneous local 
potential V(r) and considering a limited Fock space containing qq pairs only. We write the Salpeter equation 
for a meson comprising a quark and an antiquark in the meson rest frame as follows (h = c = 1): 

[sfm\ - V 2 + yjm\ - V 2 + V(r)j ip{v) = Mip{v) . (1) 

For central potential one can search for energy cigenfunctions with definite angular momentum I thus writing 

V>(r) = Y tm (t)Mr) • (2) 
The radial wavefuction 4>i(r) satisfies the equation 



dr'r 12 J dk \Jk 2 +m\ + ^ k 2 + mU k 2 j e (kr) j t {kr')4> t {r') = [M - V(r)]^(r) . (3) 



where ji(x) arc spherical Besscl functions. For I — jo{x) = smx/x and the Salpeter equation reduces to 
-J dr' J dk \ J~k 2 + ra\ + \j k 2 + rn^j sin (kr) sin (kr') u Q (r') = [M — V(r)]u (r) , (4) 

where 

u (r) = r<p (r) . (5) 

Let us discuss the potential energy. The potential V(r) we adopt comprises two parts: 

V(r) = V Ads (r)+V spm (r) . (6) 

VAd.s( r ) describes the central term of the potential, ie it contains its linearly confining part at large distances 
and the short distance behavior predicted by perturbative QCD. For VAds{ r ) we use a slightly modified 
version of the potential obtained by Andreev and Zakharov [21] in the framework of the gauge/string duality 



approach 2J, [25( . V sp i n (r) is the hyperfine term to be discussed below. 



A. Static quark potential 



In the gauge/string duality approach [2J, [25( the expectation value of the Wilson loop is given by 

(W(C))~e~ s , (7) 

where S is an area of a string world-sheet bounded by a curve C at the boundary of the AdS space (for 
references to the original papers see also 26( ) . To compute the potential the authors start with the Nambu- 
Goto action 



S = ^7 f d 2 i v/dct G nm d a X^d p X^ 
with the following background metric in D = 5 



ds A 



G nm dX n dX m = R 2 — (dx'dx 1 + dz 



l car- 
n = exp 



(8) 



(9) 
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where i = 0, . . . , 3. Choosing in © = t and £ 2 



Eq. ([5]) becomes 



S=^-T 



lx—\ll+ ( ^ 

eta; 



(10) 



where g = — -. From this action one can obtain an equation of motion for the variable z, from which an 



expression for the interquark distance r is derived as follows 2l|: 



•(A) = 2^5 £ dvv 2 exp* 1 -"'" 2 (l-^exp^ 1 -"^ 



(11) 



This expression depends on the constant c; A is a parameter in the range ]0, 2[. In terms of the original 
parameters it is given by A = cz 2 , with zq the value of the fifth coordinate z in x = 0. 

The potential is obtained by computing the energy of the configuration; first one changes the integration 
variable from x to z in (|10[) ; the resulting integral is divergent at z = and has to be regularized. The finite 
part gives 



E(X) 



9 c 

^Va 



-l 



dvv 



cxp 



Xv 2 /2 



1 — v exp 



\(l-v 2 ) 



-1/2 



1 



(12) 



This is the potential given in 21(. The dependence E(r) is obtained by elimination of the parameter A 
between (fTTj) and (p~2|) : E(r) depends on two parameters g and c. 

Our potential VAds( r ) ' m given as follows: 



V AdS (r) = 



E(r) + V 
E(r m ) + V 



r > r n 
r < r„ 



(13) 



This phenomenological procedure, which has allowed to get an expression for the QCD potential with 
the expected properties, is known as the bottom-up approach of AdS/QCD correspondence. According 
to it, starting from QCD, one constructs its dual theory in such a way to reproduce some properties of 
QCD. Besides the QCD potential, many results have been obtained so far, such as the determination of the 
numerical values of some observables |27| , in good agreement with experimental data, or the linearity of the 
Regge trajectories [28|, |29| . 



The expression (TT2")) derives from the metric (0; we have also tried to use a slight modification of the 
dilaton term in ((HJ), specifically, changing the power of z at the exponent. We find that, for any power of z, 
the potential has the same behavior, ie Coulomb-like for small r and linearly rising for large r. As far as the 
QCD potential is concerned they are all possible candidates for the AdS metric of the QCD dual space. 

We have introduced a constant term Vq in the potential. This is allowed because the more general string 
action is obtained by adding to the Nabu-Goto action a term 



5S = C R J d 2 £ R v/det G nm d a X n dpX m , 



(14) 



where R is the scalar curvature. One can easily prove that adding SS to S corresponds to add a constant 
term to the interquark potential. 



Another modification appearing in (|13p is the introduction of a cutoff at short distances r < r m [Ti 
The reason lies in the use of the relativistic kinematics embodied in Eq. (TT]). As a matter of fact, the 
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potential E(X) in (|12p diverges as A — > 0. It is a Coulombic divergence occurring at r — > 0, because the small 
distance region corresponds to A — > 0. This divergence is harmless if one uses the Schrodingcr equation 
for the wavefunction, but in the case of the Salpeter equation |T]) it produces an unphysical logarithmic 
divergence for the s— wave wavefunction for r — ► 0. To cure this unphysical singularity arising from the 
static approximation one assumes a constant potential for r smaller than 



1 

M 



(15) 



The proportionality constant can be fixed using quark duality arguments [3fJ , which shows that one expects 
different results for the equal and unequal mass cases. Thus we use 

r m = T7 (mi = m 2 ) , 

k' 

r m = J{ ( ffl l^ m l) ■ (16) 

This concludes the discussion of the central part of the potential. 

B. Spin term 

Let us now discuss the spin term. In the approximation of one-gluon exchange one knows that the spin 
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we use 



term is enhanced at short distances and is proportional to inverse quark masses. Following 

Sir) 

V spin (r) = A— ^-Si-S a (17) 
m-\m,i 

where 5 is a function enhanced at small distances. In the constituent quark model the constant A is 
proportional to a s , which is a running coupling constant, so we have introduced two different parameters 
for mesons containing at least one charm and one bottom quark (A c and Af>, respectively). We adopt the 
smeared delta function used in ref. 



31|: 



6(r)= 1^) e-" V . (18) 



The potential depends therefore on 8 parameters: c, g, Vo, k, k', A Cl At,, o. Moreover we have to fix the 
constituent quark masses: m u = md = m q , m s , m c , to;,. 



C. Numerical solutions 

To solve Eq. (fl]) we use the Multhopp method [3^], which, as shown in [l3|, is particularly useful for 
equations of the form considered here. By this method one transforms the integral equation in a set of linear 
equations introducing N parameters 9k, called Multhopp's angles. The set of equations is as follows (we 
refer the reader to 



17] for further details): 



N 



Y, B km ij(6 k ) = Myj(8 m ) (19) 
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where ip(Ok) — uq(~ cot 9 k ) , 



1 k7T 

k ~ JV+1 



with 



I jk = lim < / 

2 r (fflsmjO 



Bkr. 



d6 saijO 



(k = l,---N) and 
2 N 

3=1 



(20) 



sin 2 0(cot 6* — cot 



d9 sin j# 



2 2 
2 - - sin 9 k smj9 k 



] k) 2 + i k +e sin 2 0(cot0- cot 6 k ) 

^miA'i(mi | cot 9 — cot + (mi 7712) 



' TT J sin 2 6> L(cot6»- cot0 fe ) 2 2| cotfl - cot(9 fc | V"""^""" 

Notice that the writing in (|2ip makes explicit the prescription to avoid th< 
for the massive and the massless case; K\(x) is a modified Bcssel function. 



.(21) 



avoid the divergence in the integral both 



III. MESON SPECTRUM 



The strategy we follow to evaluate the parameters of the model is to estimate them by using information 
from the meson spectra. Subsequently, in the next sections we use this information to fit diquark masses 
and, from these data, tetraquark masses. In principle one might follow a different strategy, using all available 
experimental data to get a best fit of the parameters. We have not followed this approach because to get 
the spectrum of diquarks and tetraquarks wo make further assumptions, e.g. the use of the quark-quark 
potential as obtained by one-gluon-exchange approximation. 

Therefore, after having fixed the parameters from the mesonic spectrum, we will use below experimental 
data for X and Y states to test the hypothesis of diquark-antidiquark structure of these exotica. 

In order to get the numerical values for the parameters (c, g, Vq, k, k' , A c , Ab, a, m u = = m q , m s , 
m c , 777,5) we use the Salpeter equation for mesons with both J p = CP and 1~. We consider only mesons 
containing at least one heavy (c, b) quark, because we expect that the approximation of static potential 
works better for these states. The set of available data comprises about 20 masses. The results of the best 
fit arc reported in Table U for heavy mesons. They correspond to the following set of parameters: 



c = 0.3 GeV 2 , g = 2.75 , V Q = -0.49 GeV 

A c = 7.92 GeV 3 , A b = 3.09 GeV 3 

k = 1.48 , k' = 2.15, a = 1.21 GeV 2 

?77. 9 = 0.302 GeV , m s = 0.454 GeV , m c = 1.733 GeV , m b = 5.139 GeV . (22) 



In order to test the limits of the model we also compute the masses of states containing only light (it, d, s) 
quarks. The results for their spin averaged masses are shown in Table HT1 For the lightest mesons we find a 
large deviation, while such a discrepancy is somehow reduced in case of ss and for higher radial excitations. 
This is reasonable, since a constituent quark model with instantaneous interaction is not able to describe 
the chiral dynamics of light states. The better accuracy of the ss system allows us to fix the spin constant 
A s in eq. (TP?)) from this channel, obtaining A s =11.3 GeV 3 , which will be used in the calculation of the 
masses of the lightest diquarks. With this value, we obtain for ip a mass m=1.011 GeV (the experimental 
value is 1.019 GeV) and for ip' a mass m=1.663 GeV (the experimental value is 1.680 GeV). 

It might be useful at this stage to study the effect of the relativistic kinematics on the equation of state. 
To this end we have used the same potential, with the values of parameters indicated by Eq. (|22[) with two 



7 



TABLE I: Mass spectra for heavy mesons; q — u, d. A star (*) means that this state needs confirmation. Units are 
GeV. 



Flavor 


Level 




J — 







J 


1 






Particle 


Th. mass 


Exp. mass [33] 


Particle 


I K mass 


Exp. mass [33] 


cq 


IS 


D 


1.862 


1.867 


D* 


2.027 


2.008 




2S 




3.393 






2.598 


2.622 




35 




2.837 






2.987 




cs 


IS 


D a 


1.973 


1.968 


D* 


2.111 


2.112 




2S 




2.524 






2.670 






3S 




2.958 






3.064 




cc 


IS 


Vc 


2.990 


2.980 


J/ip 


3.125 


3.097 




2S 




3.591 


3.637 




3.655 


3.686 




3S 




3.994 






4.047 


4.039 


bq 


IS 


B 


5.198 


5.279 


B* 


5.288 


5.325 




2S 




5.757 






5.819 






SS 




6.176 






6.220 




sb 


IS 


B a 


5.301 


5.366 


B* 


5.364 


5.412 




2S 




5.856 






5.896 






3S 




6.266 






6.296 




be 


IS 


B c 


6.310 


6.286 


B* c 


6.338 


6.420 




2S 




6.869 






6.879 






3S 




7.221 






7.228 




bb 


IS 


Vb 


9.387 


9.300* 


T 


9.405 


9.460 




2S 




10.036 






10.040 


10.023 




3S 




10.369 






10.371 


10.355 




4S 




10.619 






10.620 


10.579 



TABLE II: Mass spectra for spin averaged masses of light mesons; q = u, d. Units are GeV. 



Flavor 


Level 


Th. mass 


Exp. mass [33] 


qq 


IS 


0.792 


0.616 




2S 


1.386 


1.424 


qs 


IS 


0.932 


0.794 




2S 


1.501 




ss 


IS 


0.981 


0.912 




2S 


1.571 


w 1.653 



different equations: the Salpeter equation |T]) and the Schrodinger equation. The results of the two equations 
for mesons with J p = 1~ are reported in Table HTT1 The comparison between the two computed spectra 
and the experimental one shows that, as expected, the results obtained by the Salpeter equation are more 
accurate then the ones obtained by the Schrodinger equation. We conclude that the advantage of using the 
Salpeter equation is particularly significant for the charmed states, since this equation takes into account a 
relevant source of corrections, i.e. those due to the relativistic kinematics. 

A final comment is for the use of the AdS/QCD inspired potential. As we stressed in the introduction, many 
potentials have been used in the literature to fit meson spectra. If they have to reproduce the constraints 
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TABLE III: Comparison of spectra for mesons with J p = 1 computed by the Salpeter equation in Eq. (jTJ and the 
Schrodinger equation. The same potential V(r) is used in both cases. 



Flavor 


Level 


Salpeter 


Schrodinger 


cq 


IS 


2.027 


2.154 




2S 


2.598 


2.877 


cs 


IS 


2.111 


2.182 




2S 


2.670 


2.843 


cc 


IS 


3.125 


3.133 




2S 


3.655 


3.695 


bq 


IS 


5.288 


5.494 




2S 


5.819 


6.204 


sb 


IS 


5.364 


5.507 




2S 


5.896 


6.154 


be 


IS 


6.338 


6.550 




2S 


6.879 


6.922 


bb 


IS 


9.405 


9.774 




2S 


10.040 


10.055 



of QCD they must rise linearly at large distances and have a Coulomb-like behavior at small r. We wish 
to compare the potential we have used, i.e. VacLS = E(r) + Vq in Eq. (fl3|) . with a typical QCD-inspircd 
potential, i.e. the Richardson's potential [191 ]: 



V R (r) 



Sir 



A Ar - 



/(Ar) 



33 -2n/ V Ar 
where A is a parameter, nt = 3 is the number of flavors and 

sin(qt) ( 



Vi 



4 f 

/*=-/ d 1- 

k Jo 1 



ln(l + q 2 ) q 2 



(23) 



(24) 




FIG. 1: Comparison between the AdS/QCD inspired potential VAds(r) (solid line) and the Richardson's potential 
[3 Vr(t) (dotted line). For the values of the parameters see the text. 

To allow for a comparison, similarly to what we have done for the AdS potential we have added a constant 
term V\ to the original Richardson's potential. We fix the value of A and V\ imposing that the linearly 
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confining parts of VacLS and Vr coincide at large r. For Va<is we use the fit obtained by the meson spectrum; 
this gives for A and Vi the values A = 0.44 GcV and V\ ~ - 0.6 GcV. 

We are now able to compare the two potentials and the result is in Fig. [1] It shows that, even though the 
two potentials almost coincide asymptotically, they differ significantly in the intermediate region, which is 
the region of interest from a phenomenological point of view. 

IV. DIQUARK MASSES 

In view of the applications to be discussed in the subsequent sections we wish to derive a set of parameters 
that we call effective diquark masses or simply diquark masses. 

In the one-gluon exchange approximation the color antisymmetric channel corresponding to antitriplet in 
the decomposition 

3® 3 = 3© 6 (25) 

is attractive. Moreover, in the same approximation the attractive potential between the two quarks QQ is 
half the one between a quark Q and an antiquark Q. We therefore adopt, as usual, the following value for 
the potential between the two quarks comprising the diquark: 

VQQ = \v Q Q ■ (26) 

However it must be kept in mind that this approximation holds only for the QCD perturbative interaction 
and might be modified in the linearly confining part. Therefore we expect it works better for heavy diquarks, 
as, in this case, the interaction explores smaller distances than for light diquarks. 

To derive the diquark effective masses we use the Salpeter equation with the same set of parameters used 
to fit the meson spectrum . Also in this case, as for mesons, one has that the spin term works differently 
for 5 — 1 and 5 = states. It adds a positive contribution for 5 = 1 and subtracts the same contribution 
multiplied by a factor of — 3 for 5 = 0. 

TABLE IV: Diquark masses. The masses of the present paper are obtained by the Salpeter equation (see text). The 
model in [jil. lia] uses a quasipotential of the Schrodinger type [^^.{QQ} (resp. [QQ]) means a spin 1 (resp. 5 = 0) 
diquark QQ. Units are GeV. 



State 


Mass (this paper) 


Mass [34, 351 


Mass 


Mass (this paper) 


Mass [34, 35] 


{qs} 


0.980 


1.069 


[qs] 


0.979 


0.948 


{ss} 


1.096 


1.203 








{cq} 


2.168 


2.036 


[cq] 


2.120 


1.973 


{cs} 


2.276 


2.158 


[cs] 


2.235 


2.091 


{cc} 


3.414 


3.226 








{bq} 


5.526 


5.381 


[bq] 


5.513 


5.359 


{bs} 


5.630 


5.482 


[bs] 


5.619 


5.462 


{be} 


6.741 


6.526 


[be] 


6.735 


6.519 


{bb} 


10.018 


9.778 









The results are reported in Table IIV| diquarks with 5 = are denoted as [QQ] , those with 5 = 1 as 
{QQ}. We have also reported the results of a different fit found by the authors of Refs. 3J,|35| who use a 
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Maiani and 



quasipotcntial of the Schrodinger type [36[ . We also note that in a series of papers [12 
collaborators have put forward an interpretation of the new X and Y states, as well as of light scalars such 
as cr(480), k(800), /q(980), etc., as four-quark states comprising a diquark and an antidiquark (see below). 



They use a quark constituent model similar to the well known Dc Rujula-Georgi-Glashow model [l6|. Also 



3, HI 



14. 



15| an effective diquark mass is used, however its meaning is different from ours, because 



the constituent quark mass used in these papers takes into account both the average kinetic energy and the 
potential energy, differently from our case, where a wave equation is considered. Therefore no numerical 
comparison is possible among our values for diquark masses and those found in 12, lil. [bj. fll] . 



V. TETRAQUARK SPECTRUM 



In this section we wish to discuss the possibility that two diquarks, more precisely a diquark and an 
antidiquark, combine to produce a tetraquark state. Such an interpretation was advanced long ago in order 
to give an interpretation of light scalars a (980) and /o(980) in terms of constituent quarks [37|. The recent 
discovery of new states, with both hidden and open charm, has raised a new interest for this model, though 
the interpretation of the new states is controversial, see for reviews and jij. 

Let us start with the state A(3872) @, @, 0, B 0, E3- Thc average mass of this state is 3871.2 ± 0.4 
MeV and its quantum numbers should be J PC = 1 ++ . The assignment C = +1 follows from the fact that 
the decay X — > -fj/'fy is observed. Moreover the decay X — > n + ir~ J/^> is also observed; the part of the 2ir 
invariant mass spectrum that can be ascribed to a p° decay is consistent with S— wave decay of the X state. 
From this the assignment P = +1 follows. Finally the angular distribution in this channel is incompatible 
with J = and therefore the only remaining possibility is J = 1 or J = 2. If the peak in the D°D°ir° decay 
channel at 3875.4 MeV (at 2a from the mass of AT(3872)) is interpreted as due to this state, then the J = 2 
should be excluded, which leaves us with J = 1 only. 

Several interpretations have been proposed for this state. Since the mass of this state almost coincides with 
the sum of the masses of the D° and D*° states, a natural explanation is that of a molecular state comprising 

IE HE 4(i 41, 42, 43j| . The interpretation as a x' c i state 44 is unlikely because 



the two charmed mesons 



of the small value of thc ratio 



B(X 7 J/*) 



B(X 



-J/*) 



and for the value of the mass. Also the interpretation 



as a ccg hybrid [45j is difficult due to the fact that lattice data predict larger masses for these states (for a 
discussion and for a comprehensive list of other models see 0]). 

As mentioned above, another possibility is that this state comprises four quarks. They might be a four 
quark cluster without internal structure 46| or a bound state of two diquarks, as first discussed in [l3| and 
more recently in 35[ . We shall return to this interpretation below. 



Thc Belle collaboration observes two bumps, at 3940 MeV. They do not necessarily correspond to two 
different states. Thc state called F(3940) is observed in the decay mode B — > KloJ/^ (ll|. Its reported 
mass is M = 3943 ± 11 ± 13 MeV. Its interpretation as a charmonium cc state is possible, but should be 
corroborated by the observation of the decay mode Y — > D^D^*\ which has not yet been seen. Also the 
interpretation as a cc-gluon has been proposed, though the predicted mass of such a state, around 4.3-4.5 
GeV from lattice QCD computations 



47 



is significantly larger than the measured value. Also in this 
case a four quark interpretation is possible, in particular the state might be the 2 ++ state predicted by the 
diquark-antidiquark scheme [IE]- The other state A(3940) needs confirmation; it is observed 49| in double 
charm production: e + e~ — > J/tpX — » J/ipDD*. Its possible interpretations are the states x' c 01 Vc 



We wish to test the interpretation of these states as diquark-antidiquark bound states, in thc same spirit 
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ofRefs. [II El, El and [34(. Note however that our approach is different from that followed in these 
papers. As a matter of fact we apply a wave equation for the tetraquark states, comprising the diquarks 
(Q1Q2) and (Q3Q4), as follows 



'12 



V 2 



'34 



V 2 + V(R) MK) = MtVt(R) 



(27) 



Here M t and ipt(R) are the tetraquark mass and wavefunction respectively, m.y is the effective diquark 
mass computed above; R is the distance between the centers of the two diquarks. We take into account the 
structure of the diquarks by defining a smeared potential as follows: 



R + ri — r 2 



with 



N 



dri / dr 2 |V>i2(ri)| 2 |<Mr2)P 



(28) 



(29) 



In these equations "012(1*1) (resp. ip3i( r 2)) is the wave function of the diquark whose content is Q1Q2 (rcsp. 
Q3Q4); for V(r) we use ©, and the masses appearing in the spin term of Eq. (|17p arc diquark masses. In 
this paper we only consider diquark with internal orbital quantum number £ = 0; therefore tpij is related to 

Uo(r) defined above as follows: "/>ij(r) = U ° ■ 



0.25 
0.2 

0.15 
0.1 

0.05 



FIG. 2: The squared diquark wavefunction (in GeV 3 ) versus r (in GeV 1 ). Data correspond to the diquark [eg]. 

Since |i/'i2(r)| 2 is strongly peaked at r ~ (see Fig. [2]) we cut-off the integrals in Eqns. (|28|) and (|29|) at 
t~i,2 < -Ro where i?o is the peak value of Mot(r), the wavefunction of the tetraquark. This procedure ensures 
that the two diquarks are on average inside the tetraquark's bag. In Fig. [3] is represented Mot(r) for the 
tetraquark [eg] {eg}. 

We finally note that we use for ^12 (r) the result obtained from the diquark wave equation. This is only 
approximately correct because that equation provides the diquark wavefunction in the diquark rest frame, 
whereas (|27[) holds in the tetraquark rest frame. However for diquarks comprising heavy quarks (c, b) the 
average diquark velocity is small (we estimate j3 ~ 0.15 for diquarks with open charm and j3 ~ 0.06 for 
diquarks with open bottom). Therefore we can neglect the distortion induced by the Lorentz boost on the 
wavefunction. The effect of convoluting the potential with the wavefunctions in (|28p is shown in Fig. 3J 

In Tables [V] and I VII we present our predictions for the four-quark states with hidden charm and hidden 
bottom respectively. Notice that we have adopted the same values of the free parameters used in the 
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FIG. 3: The tetraquark wavefunction uot(r) (in GeV 1//2 ) versus r (in GeV 1 ). Data correspond to the tetraquark 
[cq]{cq}. 

V 

1 r 




FIG. 4: The potential between static diquarks (dashed line) and its modification according Eq. (|28|) (solid line). 
Units are. GeV (V) and GeV -1 (r). Data refer to the [cg]{cg} potential. 



previous sections, so basically our calculation is parameter-free. It is remarkable that in the two cases 
when a comparison with experiment is possible, our results agree, within the errors, with the data. A 
further comment is the following. A peculiar feature of the diquark-antidiquark scheme for the X state 
is the prediction of four different states with mass differences of a few McV. Two of them arc neutral: 
X u = [cu][cu], Xd = [cd][cd], and two charged: X + = [cu][cd\, X~ = [cd][ct2]. It is possible that the two 
states seen through the decay modes J/ipir + ir~ and D Q D°ir° are really different and they would correspond 
to the neutral states X u , Xd since the mass difference between the two bumps at 3871 and 3876 MeV is of 
the same order of magnitude of the mass difference between the u and d quarks. According to [lil this is a 
piece of evidence in favor of the four quark interpretation of these states. We have not included any mass 
difference between u and d quarks, but also in our scheme one would expect a mass difference of a few McV 
between the two neutral states. We note that in lj| a strategy to find the two missing charged partners is 
delineated. We refer to this paper for further details. 

By our model we can compute also radial excitations of tetraquarks. For example the first radial excitations 
of the two X states with 1 H have mass m = 4.421 GeV and m = 4.418 GeV respectively. In Ref. [5(j] the 
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TABLE V: Four-quark states with hidden charm interpreted as bound states comprising a diquark (cq) and an 
antidiquark (cq). {QQ} (resp. [QQ]) means a spin 1 (resp. spin 0) diquark QQ. The masses of the present paper 
are obtained by the Salpeter equation (see text). The model in [3, EHl uses a quasipotential of the Schrodinger type 
0. Ref. uses a constituent quark model (' means that the experimental value is used as an input in this case). 
Units are GeV. 





Flavor content 


Mass (this paper) 


Mass [34. 351 


Mass [131 


Exp. State 


Exp Mass 


0++ 


[cq][cq] 


3.857 


3.812 


3.723 






1++ 


({cq]{cq} + [cq]{cq})/V2 


3.899 


3.871 


3.872+ 


X(3872) 


3.8712 ± 0.0004 [5, 6, £ .8, 9, 10} 


1+- 


({cq]{cq} - [cq]{cq})/V2 


3.899 


3.871 


3.754 






0++ 


{cq}{cq} 


3.729 


3.852 


3.832 






1+- 


{cq}{cq} 


3.833 


3.890 


3.882 






2++ 


{cq}{cq} 


3.988 


3.968 


3.952 


V(3940) 


3.943 ±0.011 ±0.013 [11] 



TABLE VI: Four-quark states with hidden bottom interpreted as bound states comprising a diquark (bq) and an 
antidiquark (bq). {QQ} (resp. [QQ]) means a spin 1 (resp. S — 0) diquark QQ. The masses of the present paper 
are obtained by the Salpeter equation (see text). The model in [3, EH] uses a quasipotential of the Schrodinger type 
Units are GeV. 



J PC 


Flavor content 


Mass (this paper) 


Mass [34, .3J 


0++ 


[bq][bq] 


10.260 


10.471 


l +± 


([bq]{bq} ± [bq]{bq})/V2 


10.284 


10.492 


0++ 


{bq}{bq} 


10.264 


10.473 


1+- 


{bq}{bq} 


10.275 


10.484 


2++ 


{bq}{bq} 


10.296 


10.534 



state Z(4433), recently observed by the Belle Collaboration [5l| through the decay Z(4433) — > ijj(2S)Tr , is 
interpreted as the first radial excitation of one of these states. Although the difference between theory and 
experiment in this case is larger than, say, for X(3872) or F(3940), this interpretation is compatible with 
our results because of the theoretical errors of the present model. Another exotic state is F(4260), found by 



the BaBar Collaboration [52J and confirmed by CLEO [53j and BELLE [54[. It is interpreted in [14| as an 
orbital excitation of a tetraquark state, an interpretation we are neither able to confirm nor to disprove, as 
we have limited the analysis to the I — states. 

Let us now comment on the light tetraquarks. We do not include them in the Tables because for them the 
assumptions of our model and therefore its results are less reliable; moreover, as discussed in connections 
with Eqns. I|28p. (|2T)|) . the distortion of the diquark wavefunction due to the relativistic motion is larger for 
light diquarks and its neglect generates a greater error. 

Let us finally comment on the possible existence of tetraquarks comprising a heavy diquark and a light 



diquark. We present in Table IVIII our predictions for tetraquarks with open charm and strangeness and 
compare them with the prediction of the constituent quark model of Ref. [l3j]. In fill th e state + is 
associated with the particle £> s (2317) 0, 1 + with D s (2457) Q and 2+ with X(2632) Q. Our results 
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TABLE VII: Comparison between the results of the present model and those of Ref. [l3|] for tetraquarks with open 
charm and strangeness. Units are GeV. 



J p 


Flavor content 


Th. mass (this paper) 


Th. mass (model [131*1 


0+ 


[cq][qs] 


2.840 


2.371 


0+ 


{cq}{qs} 


2.503 


2.424 


1+ 


{cq}[qs] 


2.880 


2.410 


1+ 


{cq}{qs} 


2.748 


2.462 


1+ 


[cq}{qs} 


2.841 


2.571 


2+ 


{cq}{qs} 


2.983 


2.648 



are significantly different from those of [13|. Again, this might be due to the limitations of one or both the 
constituent quark models. In any event, on the basis of the results in Section HTT1 we do not expect theoretical 
errors larger than a few hundred MeV for the results of the present model in Table IVIIl so that we do not 
support the interpretation of the states Z? s (2317), £> s (2457) and A(2632) as tetraquark charmed states with 
open strangeness. In [5lj ] these states are interpreted as a mixture of P-wave quark-antiquark states and 
four-quark components. 



VI. CONCLUSIONS 



We have developed an application and, at the same time, a test for the QCD potential found by means of 
the AdS/QCD correspondence. We have put it in a semirelativistic wave equation and fitted meson spectra. 
Our result is that this model, with the AdS/QCD potential plus a contribution from spin interaction, can 
reproduce the experimental data except for the lighter states (n, K). This agreement has motivated us 
to make some predictions on the masses of tetraquarks, considering them as bound states of a diquark 
and an antiquark. Our conclusions are that some tetraquark states with appropriate flavor content can be 
identified with the particles A(3872) and y(3940). On the other hand the present model does not favor the 
interpretation of some charmed positive parity particles with strangeness as tetraquark states. 
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